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Abstract 

The optical implementation of complex network is considered for 
laser cavity formed by active mirrors deposited on the surface of plane- 
parallel slice, ellipsoidal microball , whose eccentricity is close to 0.5 
and equivalent optical circuit, containing amplifying and absorbing 
media in Fourier-conjugated planes. In this types of microlaser's the 
self-imaging cavity provides efficient diffractive interconnections be- 
tween different points of the mirrors. The interconnections could be 
local, as in plane-parallel geometry, global as in confocal cavity or 
inhomogeneous one with definite subset of nodes, having much more 
interconnections than their neighbours. We show how the spatiotem- 
poral dynamics of this nonlinear resonator is controlled by the geom- 
etry of gain-losses distribution making possible the analog modeling 
of the several generic nonlinear evolution equations. 

PACS number(s:) 05.45-a., 42.60.-v, 45.70.Qj 



1 Introduction. 

The close analogy between spatiotemporal patterns in optics and hydro- 
dynamics is being the subject of relentless interest within recent decades 
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[H El El S] • Inspite of a substantial difference between Maxwell-Bloch and 
Navier - Stockes equations, the definite similarity of dynamical behaviour 
have been already established O |6] . It is well known that both set of equa- 
tions are reduced with the some degree of accuracy to the Ginzburg-Landau 
(GLE) [3] or to the Swift-Hohenberg (SHE) |6] equations. The last equation 
contains additional fourth-order dispersion terms, which arrests collapse, in- 
herent to 2D GLE and stabilizes the spatial solitons [B]. The competition 
between localized structures [71 [8] , stripes, lattices and chaos is shown to be 
highly sensitive also to boundary conditions, which could be controlled by the 
proper choice of the Fourier spatial filter [9] placed inside the cavity. In con- 
trast to the case of spatial filtering by finite gain linewidth [B], the boundary 
conditions itself, namely the spatial layout of the gain elements, diaphragms 
and mirrors provides the possibility to control the spatiotemporal dynamics 
[H H]. In this context the situation in laser cavity is similar to the complex 
networks physics, where spatial interconnections define both dynamical and 
statistical properties of network [TT] . 

It is shown in a set of interesting cases that conventional set of Maxwell- 
Bloch equations for class-A laser could be reduced by split-step Green func- 
tion technigue to the integral equation including the boundary conditions in 
explicit form [H [10] : 



/ K{r-r')f{En{r'))(fr' (1) 

-oo J —oo 

This equation has a form of the product of two operators in the form of 
convolution: 

E^+,{r)=kf{E^{r)) (2) 

where is amplitude of electromagnetic field, K is linear nonlocal operator 
(propagator of the cavity), / is nonlinear local operator obtained under thin 
slice approximation [^. For example, the kernel K for the perfectly plane- 
parallel Fabry- Perot cavity of microchip laser has the form p!2l [13] : 



—* ik — * 
K{f—r') = exp{ik{f — r'Y / 2z) (3) 

the transfer function f[En{r),G{r)\ explicitly contains the gain distribution 
G{f) for this microchip laser cavity configuration (fig.l) [TOl [T3l [TH [T5] : 

On the other hand, for the perfectly confocal cavity (fig. 2), the kernel K 
is Fourier transform of diaphragm transmission. For example, when confocal 



2 



R ^ 0.999 



OutputMirror 
i?^ 0.97 



X 



d ~ 500/im 
Nd : YAG 
slice 



osc 




Intensity 
2D 
vortex 
array 



Figure 1: High-Q microchip solid-state laser cavity p3] . The length z = 
is close to beam diameter D. The transversal r = (x, y) intensity distribution 
in 2D vortex array inline. 
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Figure 2: a)Quasiconfocal cavity formed by active mirrors on the surface of 
the ellipsoidal microball[IO] . The focal length of active mirrors F is a half 
of raduis of curvature Rc ■ b) Quasiconfocal cavity with two lenses and two 
plane-parallel nonlinear elements: (^(r") and f{En{r)) [16]. 
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cavity is formed by two mirrors, located in Fourier - conjugated planes, one 
being active, the kernel K becomes the Fourier transform of gain distribution 
G HIIIO]: 

K{f-r')= / exp{ik{f-r')r"/2z)G{r")d^r" (4) 

J —oo J —oo 

The confocal cavity could be formed by two concave mirrors with the common 
focus [To], four mirrors in 8-F ring laser geometry ^ or inside the ellipsoidal 
microball (fig.2a). The equivalent optical circuit formed by two confocal 
lenses and two planar, gain and absorber layers, had been used for excitation 
and detection of the spatial solitons (fig. 2b) In the plane-parallel geom- 
etry (fig.l) the interconnections occur mainly in between adjacent points of 
the mirror (neighbour-neibourhood interaction), while in confocal geometry 
(fig.2) the interconnections are global: each point interacts with the others 
through gain layer, located in "far field" . It is clear enough that the in- 
termediate situation in between plane - parallel (fig.l) and confocal (fig.2) 
configurations, is described by equation ([T]) as well, provided the kernel K 
contains terms corresponding to both difraction in free space and absorbance 
at the diafragm ("diffusion"). This intermediate situation with possibility of 
creating the selected "nodes" having much more interconnections than neigh- 
bours is interested from the point of view of complex networks theory [TT] . 

We will see that initial master equation ([T]) having the form of nonlinear 
nonlocal map, the nonlocality being induced by the explicit inclusion of the 
boundary conditions [1] , is a quite general form of intracavity laser dynamics. 
Originally this model had been used for modeling of travelling wave tube 
microwave oscillator wih delayed feedback and some other relevant systems 
[T7j . the scaling relations for a weak chaos had been established. Afterwards 
the same model, called "continuous coupled map" (or "continuous coupled 
map lattice"), had also been used for description of the pattern formation 
in vibrating granular layer |18]. Our goal here is to show that nonlinear 
dynamics described by primary master equation in the form of nonlocal map 
([T]) contains the different types of nonlinear behaviour, intrinsic to a wide set 
of several nonlinear evolution equations. The dynamics inherent to Burgers, 
nonlinear Shrodinger, Ginzburg-Landau or Swift-Hohenberg equations will 
be released by the proper choice of gain-losses distribution. 
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2 High gain limit. 



In the high gain hmit, when the field changes in nonhnear shce are assumed 
to be arbitrarily large compared to diffraction effects, we have : 



f{En{r)} + ai — — h a2 



dE dx 



dE dx"^ dE"^ dx 



+ 



+ "3 



dfd^En , d^f dEnd'^En , dp fdE„ 



where 



a 



dE dx^ 
1 



+ 3 



dE"^ dx dx"^ dE^ \ dx 



{x - xTK{x - x')dx' 



+ ...(5) 



(6) 



are the moments of kernel K, the one transverse dimension x is considered 
here and further for simplification of expressions. Now we ought to take into 
account the components of the gain function / having different orders of 
magnitude : 



f{En{f))=GdEn l-Ps\En 



(7) 



Obviously, this form of the gain function corresponds to the two-level medium 
with weak saturation. When 02, as = we have the wave equation with 
nonlinear source f{En{r)): 



E^+,{f) = f{E^{f)) + 



dE 



dEn 
dx 



(8) 



The transition from "discrete time" n with step 6t 



2c 



is the speed of 



light ) to continuous time t is straightforward and we can get in the so-called 
"mean-field approximation" O [6l [19] the equation for the order parameter: 



dE 



dE 
dx 



(9) 



where Tc is cavity lifetime. Note the similar equation appears also in acoustics 
[20j . When as = the Burgers-like equation with nonlinear source / follows: 



ai 



dl 

dE 



dK 
dx 



+ 



' ~ Of 
""'dE 



d^E„ 
dx^ 



(10) 
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or for continuous time t: 
dE{f, t) 



dt 



E 



+ 



dE 



ai 



dE 
dx 



+ 



dE 



0L2 



d'^E 



(11) 



When a2 = the analog of the Korteveg-De-Vries equation (again with 
nonhnear source/ ) occurs: 



dE 



dx 



+ 



"3 



dE 



d^ 
dx^ 



or for continuous time t 
dE{f,t) 



E 



dt -f(E(r,t))- + 



~ Of 
""'dE 



dE 
dx 



+ 



' ~ Of 
""'dE 



d^E 
dx^ 



(12) 



(13) 



3 Low gain limit. 

Now turn our attention to the low gain limit: nonlinearity and dispersion are 
of the same order 6 . We are applying here the same "formal" trick as in [1] 
using the following identities: 

K{x - x') = 6{x - x') + {K{x - x') - 6{x - x')} (14) 

" V ' 

and 

/(E„(x)) = E^{x) + {/(^„(x)) - Enix)} (15) 
« ' 

In f|T4|) and f|T5l) the small terms underbraced by figure brackets, are re- 
sponsible for nonlinearity and dispersion. They are assumed to be of the 
order 6 of smallness . Hence we have another secondary master equation, by 
substitution of fflil) and (fT5|) to (fT]) : 

E„,^,{r) = f{E^{f))-E^{r) + kEr,{f^ + 

+ {kE^{x)-E„ix)}{f{En{x))-En{x)} (16) 



Because the last term in f|T6l) is of the order 6^, it could be omitted and we 
have the secondary master equation in the following form: 

En+i{^ = /(K(r)) - Er,{f) + KErXr) (17) 
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and for continuous time: 

= f{E{f, t)) - 1 + hm t)) (18) 

Applying the same series expansion to the convolution inside (1171) we have: 

oo om Tp 

E^^,{r) = /(E„(f)) - E„(r) + Y: o^m-j-^ (19) 

m=l '-'^ 

where am are again the moments ([6]) of the kernel K. Consequently for 
continuous time one may write: 

dE{f,t) E ~ 9™E , , 

4^ = /(B(f.*))--+E"™a^ (20) 

As a further approximation one can get readily the several familiar dy- 
namic equations, depending on relative magnitude of real and imaginary 
parts of the moments of kernel and relative strength of the nonlinearity. For 
example for all am = except for m = 2 and a2 - real (pure confocal cavity), 
Kolmogorov-Petrovsky-Piskounov equation follows (KPP) [11 [TO]: 

E^+,{x) = f{E^ix)) - E„{x) + (21) 

For continuous time this equation has the following form: 

dE(f,t) E d'^E 

Next, for the ai = , a2 - purely imaginary (plane-parallel Fabry-Perot cav- 
ity), and the /(i?„(x)) = ikn2dEn\E\'^ nonlinear Shrodinger equation (NLS) 
occurs [3l H] : 

d'^E 

En+i{x) = ilm{a2)-^^ - -E„ + ikn2dEn\En\^ (23) 

For continuous time we have: 

dE{r,t) .^.d^E E 



ilm{a2)^r-ir \- ickn2E\E\'' (24) 



dt dx 
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Furthermore, for all moments of kernel K am = except for m = 2 , and 
a2 - arbitrarily nonzero complex number ( Fabry-Perot cavity with curved 
mirrors) and 



f{Enix)) = GdEn 1 + {ReP)En\En\^ + tiImP)En\En\ 



(25) 



(semiconductor or solid-state gain medium) the Ginzburg-Landau equation 
(GLE) is valid [4^: 



E, 



n+ll 



Re[a2) ^ „ + ilm[a2) 



dx"^ 



EJx) + 



+GdEn 1 + (i?e/?)E„|E„|2 + i{ImP)En\E„\ 



(26) 



For continuous time we have: 



dE{r,t) 
di 



Re{a2 



^.d^E E 



dx"^ Tr 



-ilm{a2 



dx"^ 



-GE+{Re(3)E\E\'^+i{Im(3)E\E\'^ 

(27) 

In a final step of our consideration, for the only nonzero a2 and - arbitrarily 
complex numbers ( corresponding to Fabry-Perot cavity with curved mirrors) 
and 



f{En{x)) = GEn + /?^„|^„r + lEn\E^ 



(28) 



( semiconductor or solid-state gain medium in deep saturation ) the Swift- 
Hohenberg equation (SHE) tooks place [20] : 



En+i{x) = a2^-jr + «4^rT- 



dx"^ dx^ 
And at last we have for continuous time: 
dE{r,t) ^d^E ^ d^E E 



Enix) + GEn + pEn\En\^ + 7K|^n|^ (29) 



dt 



a2 



^ + ^^^-- + OE + mE\^ + ^ElEl 



(30) 



4 Conclusion. 



We have just saw that possible dynamical regimes of the electromagnetic 
field in microlaser cavity is much wider than those simulated by parabolic 
or Swift-Hohenberg equations scenarious. The boundary conditions itself 
affect the form of the master equation, hence they affect the observed laser 
dynamics. 
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It is noteworthy to keep in mind that each of the evolution equations 
derived above have its own pecuhar behaviour. Some of these equations, hke 
KPP fl2T|) or NLSE fl23l) have exact solutions in the form of kinks (switch- 
ing waves ) or solitons . As for the other equations, like GLE fl26l) or SHE 
fl29|) . their exact solutions are not known, although there exists a set of 
established analytical, semianalytical and numerical results, showing which 
dynamics is to be expected in the given range of parameters. Moreover by 
means of choosing the special resonator configuration, setting the appropri- 
ate gain-losses distribution G and nonlinearity it is possible to simulate the 
spatiotemporal dynamics of a given dynamical equation , say GLE (!26|l or 
SHE equations (1291) . The initial conditions will evolve along the phase trajec- 
tories of the evolution equation to be investigated. This provides a possibility 
of the purely optical analog modeling of partial differential equations ([H [121 
|26| [211 [23| [26l [29]) extracted above from nonlocal maps (coupled map lat- 
tices). On the other hand because this approach is going to be useful also 
in condensed matter physics, for example as a model of granular materials 
|18j : the results presented here provide the definite information concerning 
the expected dynamical regimes of the pattern formation and properties of 
this materials. 
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